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^ As a particle moves through a fluid, it may generate a laminar wake behind it. In 

the gauge-string duality, we show that such a diffusion wake is created by a heavy quark 

^ moving through a thermal plasma and that it has a universal strength when compared to 

^ the total drag force exerted on the quark by the plasma. The universality extends over all 

asymptotically anti-de Sitter supergravity constructions with arbitrary scalar matter. We 
discuss how these results relate to the linearized hydrodynamic approximation and how they 
bear on our understanding of di-hadron correlators in heavy ion collisions. 
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1 Introduction 

Consider a quark moving at constant velocity through an infinite, static, thermal medium. 
Assume that the temperature is above the confinement transition, so the medium is a plasma 
of gluons and whatever dynamical matter there may be. The quark loses energy to the 
plasma, but where does this energy go? Sufficiently far from the quark, a complete account 
of energy dissipation should be possible within the framework of linearized hydrodynamics. 
Hydrodynamics is usually applicable at length scales much larger than the mean free path 
of thermal quasi-particles within the plasma, which is to say some multiple of 1/T. (We 
set h = c = = 1 throughout.) The applicability of linearized hydrodynamics depends on 
being far enough from the quark for the perturbations to be very small. It is well understood 
(see for example [1]) that if one assumes axial symmetry of the fluid flow around the direction 
of motion of the quark, then the equations of linearized hydrodynamics partially decouple 
into a pair of equations describing compressive waves (sound) and a third equation describing 
shearing modes (the diffusion wake)]^ Figure [l] shows a typical geometry for these two modes. 

^In principle, there could be additional modes coming from conserved currents. In the scenarios we 
consider, such currents do not appear. 
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We denote the total energy of the moving quark as E. The total energy gained by the 
plasma is given by P* = —dE/dt, where in our conventions dE/dt < means that energy 
is flowing from the quark into the plasma. Assuming that sound and diffusion are the only 
modes of energy loss from the quark, the total power delivered to the plasma is 

pt =pd^ps^ (1) 

where and P'^ correspond respectively to the energy lost to sound and the energy lost to 
diffusion. 

We claim that in a fairly broad range of string theory constructions and using a sensible 
but non-unique way of separating diffusion modes from other modes, 

: pt = -1 : . (2) 

The relative minus sign means that the diffusion wake supplies energy to the quark. We 
prove this in section |2] for any theory admitting a holographic dual asymptotic to AdSs that 
can be described in terms of Einstein gravity coupled to scalars with no more than two 
derivatives in the action. In our demonstration, the quark has to be infinitely massive, and 
it is represented as a string trailing into the five-dimensional dual geometry]^ The tension of 
the string is an arbitrary function of the scalars. 
Putting ([T]) and ^ together, we learn that 

P' : P'^ : P' = l + v^ : -l-.v"" . (3) 

This means that energy flows into sound waves at a rate which is 1 + l/f^ times as fast as 
the total outflow of energy from the quark. The diffusion wake supplies energy to the quark 
in just the right amount to balance equation ([T]). In section |3] we show that we do not need 
to assume ([T]) in order to demonstrate ([s]), at least in the case of a single scalar. The result 
(|3| holds for all v between and 1, but the most interesting case is when v is greater than 
the speed of sound Cs, in which case there is a sonic boom. 

The result Q was demonstrated for A/" = 4 super- Yang-Mills with a large number of 
colors and a large 't Hooft coupling in |2], which followed upon a number of works by 
ourselves and other authors [31 IH El El El El E] ; parallel lines of development have included 



^Clearly an infinitely massive quark must also have divergent energy E. But P* = —dE/dt is finite and 
closely related to the drag force as calculated for A/" = 4 super- Yang-Mills theory in [21 H]; see also the 



discussion at the end of section 2.1 
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Figure 1: Top: A cartoon of the diffusion wake and the sonic boom. The Mach cone and 
the diffusion wake experience comparable viscosity broadening, controlled by the parameter 
Ts = Arj/STs. Bottom: The Poynting vector 5* produced by a heavy quark propagating 
through a thermal plasma of A/" = 4 super- Yang-Mills theory, from [2]. The arrows indicate 
the direction of S, and the color indicates its magnitude, with red meaning large IS"]. The 
dashed green line marks the Mach angle, and the gray line marks the curve where the far 
field asymptotics of the wake falls to half its maximal value. 
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[iniiniinillSlElIISlIISlin]. We note in particular the close relation of [12] to the drag 
force results of [31 H] . 

The backgrounds we consider can accommodate arbitrary speed of sound (limited, per- 
haps, by Cs < including a speed of sound that varies as a function of temperature. 
Thus our result ([s]) has a degree of universality which encourages us to think that it might be 
compared with data from heavy ion collisions. As we explain in section [5} a diffusion wake 
as strong as we find poses a challenge to our understanding of the experimental situation. 

2 The diffusion wake 

The action we start with is 



^5 

where is the scalar Lagrangian 



2^5 Jm 27ra' Jw 



(4) 



= -\G^''nu{$)d,<p'dA' - V{4>) ■ (5) 

G^y is the metric on the bulk spacetime M, and gap is the metric induced from it on the 
string worldsheet W . We will consider backgrounds which are asymptotically anti-de Sitter, 
with (p ^ (pB near the conformal boundary of AdS. The functions V{(f)), and g(0) are 

arbitrary, except for the requirement that V{(f)) should have a local (negative) maximum 
for a finite value of the scalars, V{(j)B) = —12/-^^, in order for the desired anti-de Sitter 
asymptotics to arise near the conformal boundary. 

When the string is not present, the backgrounds of interest take the form 



ds = air 



\2 



-hir)dt^ + dx^ + — - 
h[r) 



(6) 



This ansatz is the result of assuming translation invariance in the R^'^ directions and rota- 
tional invariance in the directions. It is assumed that /i 1 at the conformal boundary 
where also r — and a L/r. There is an event horizon at r = rj^ > where h{rii) = 0. 
The tt and ii components of the Einstein equations (where i runs over the three components 
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of x) provide us with 



2 ha a 2 h 
where primes denote d/dr. The scalars satisfy 



3 h'a' a" 1 

+ 3— = -^i:<^, (7b) 



d dLtf, dL^ 



. (7c) 



drdcj)^' d(j)' 

In addition, the rr component of the Einstein equation imposes a zero-energy constraint: 

3 h'a' a'"^ la^ , ^ <r,<j,\ /on 

The system of equations ([T]) and ([s]) has been extensively studied, and we briefly sum- 
marize aspects of the relevant literature. It was shown in flF] that a solution of the form ^ 
can exist in the case of a single scalar only if V{(j)H) < ^(0b) with (pn = 4'{fH)- Numerical 
construction of such black holes was first reported on in [TH] for a specific potential. Further 
numerical solutions were exhibited more explicitly in [T^ and later in [201 Ell 1221 [221 [21] , 
where analyses of thermodynamic quantities were also given. Following [25], some approxi- 
mate, analytical results for high temperature black holes were obtained in [261 [2Z1 128] . 



2.1 The trailing string 

In [31 [1] the drag force acting on a moving quark at constant velocity v through a plasma has 
been calculated using AdS/CFT. It is straightforward to generalize the drag force calculation 
of [31 H] to the backgrounds described above. Various generalizations exist in the literature, 
for example [291 IBOl [3TJ. The one we will describe here is quite a modest extension of 
previously published results, which allows for an arbitrary dependence on the scalar field 
(j)^ coming from a Kaluza-Klein reduction of the ten dimensional Dilaton and metric. The 
ansatz for the trailing string is 

= vt + ^(r) , (9) 

and the string action takes the form 



S=ldtdrL L = -^a\l-'^ + he. (10) 
' 2Txa' y h 
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The conserved momentum conjugate to ^ is 



The drag force exerted by the plasma on the quark is 

f, = =^ = -.,. (12) 
at 

Note that tt^ > 0, so Fi is negative, meaning in a direction opposite the quark's velocity. 
The total rate of energy loss from the quark to the plasma is 

pt dE dE dpi , . 

at ax at 

which is a positive quantity. 

2.2 Perturbations in axial gauge 

The trailing string sources the metric and the scalar fields coupled to it. Treating as a 
small parameter, we may evaluate the response of the various fields to the string to linear 
order in nf. 

We consider linearized fluctuations around a background of the form ([6]). These fluctua- 
tions will be expressed as 

Gf,^ ^ G^^ + 6G^^, 0^0 + 50, etc. (14) 

In axial gauge one sets 6Gru = for all u. We further define 

5Gmn = nla^Hjnn , (15) 

where the pre-factor KgO;^ is for convenience. 

As in [32j, we study the linearized fluctuations in Fourier space with a co-moving ansatz: 

Hmn{t, r, f) = j ^e^'=i(^^-^*)+'=^^'+'^^^'ff^„(t, r, k) . (16) 

Because of axial symmetry, the Fourier coefficients Hmn can only depend on wave-number k 
through ki and k± = \/k\^ k\. Let's fix this symmetry by setting fcs = and k2 = k± > 0. 
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The three components {Hq^, H13, and -^23) are odd under reflection through the plane 
perpendicular to 3. The trailing string is symmetric under this reflection, so it doesn't source 
the odd metric components. We may therefore set them to zero. The other components [Hqo, 
Hqi, Hq2, Hii, H12, H22, and -^33) are even, and they are all sourced non-trivially by the 
trailing string. A substantial simpliflcation of the equations of motion arises from forming 
the following linear combinations, which are essentially the same as the ones in [32] : 



A 




-Hii + 2-^Hi2 



H22 + 



k 



H- 



33 



(17) 



D 




1 

2^2 



h [Hoi 



02 



fcl 



H12 + H- 



22 



E 







E2 













\ {—\Hqq + Hii + if 22 + -^^33) 

2b (^01 + 1^Hq2 



\ 



\ {Hii + H22 + -ff33) 

\\ {Hu (2 - 6 (^)^) + H22 (-4 + 6 (^)^) + 2if33 - 12^//i2) ) 



■ (19) 



The equation of motion for A decouples from all other metric components as well as from 
the scalars. The equations of motion for D mix with one another but decouple from A and 
E as well as the scalars. The equations of motion for E couple to the scalars. The deflning 



relations (|17|)-(|19|) are equivalent to 

-El + E3) 



Hoo 
Hqi 

H02 
Hn 
H12 

H22 
H33 



2h 



2v 

2vkik_i 
k^ 

v^k^, 



^^{klDi + klE2) 



-Di + E2 



Av^klkl 2 -2kl + kl 
W"^ - -1^^' + 3^^ + 3P 



(20) 



kiki 



v^k"", . ^ ,kl 



F V k^ 
v^klk^ 



■-A + 2v' 



-Do - E, 



A;2 

^ Av^klklD2 2 ^ kl-2kl^ 
P + 3^^ + ^fc^^^ 



— + -E3 + -Ei . 
fc2 3 3 



It will prove useful to have these inverse relations ready to hand when we want to extract 
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components of the gauge theory's stress energy tensor. 

We shall restrict our attention to the vector modes since these are the only ones which 
will contribute to the diffusion wake. To see this consider the contribution of the moving 
quark to the Poynting vector 

S=(^Si S2 53) where S, = = -{6Toi) . (21) 

Recall that we employ mostly plus signature and that (T"*) = in the static background. 
We use S in preference to 6S for simplicity, keeping in mind that the non-zero Poynting 
vector owes wholly to the presence of the moving quark. Using holographic renormalization 
(see appendix |A]) we can show that 

S, = - hm ^i/o. . (22) 

It is useful to decompose the Fourier components of the Poynting vector into a longitudinal 
part and a transverse part 

S = SLk + ST, (23) 

where k = j:{ki,k±). The St component is orthogonal to k and lies in the same plane as 
k and v due to the azimuthal symmetry around the direction of motion of the quark. It is 
well understood (see for example pj) that the diffusion wake is described in terms of St- 



Referring to (20), we see that the contribution of E2 to 5* is in a direction parallel to the 



wave-number k, whereas the contribution to Di is orthogonal to k. In short, E2 controls St 



while Di controls St- Using (20) and (22) one finds 



The linearized Einstein equations in axial gauge imply the following second order equa- 
tions for the D modes, 

(25) 



In (25), the real- valued function ^(r) is determined by solving (11). Additionally there is a 
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first order constraint, 



D[ - hD'^ 



(26) 



Note that the source for the constraint equations (26) is of order 0{k ) while that of the 



second order equations (25) is of order 0{1). 



2.3 Large distance asymptotics 

We are interested in the large distance asymptotics of the energy momentum tensor. In this 
regime one would expect that linearized hydrodynamics will be a good approximation. The 
equations for vector perturbations simplify drastically when only the leading order terms in 



small k are retained: the equations of motion (25) become 

d 



dr 







d 
dr 



and the first order constraint (26) becomes 



D[ - hD'^ 



(27) 



iv'^kio^ 



(28) 



The most general solution of the equations of motion (27) consistent with the boundary 
condition Di 



at the conformal boundary of the asymptotically AdS background is 
df „ , X , r dr 



DAr) = di 



a(r)^ 



D2{r) = d2 



Q /i(f)a(r)3 



(29) 



Because /i — at the horizon, D2 has a logarithmic divergence there. Requiring that the 
metric fiuctuations are finite at the horizon implies d2 = 0. This is equivalent to taking the 
k ^ limit of the usual horizon boundary condition where outgoing modes are forbidden. 
A more detailed discussion of the horizon boundary condition can be found in section [3] 



Plugging (29) into (28) with ^2 = leads immediately to 



di 



iv'^ki 



(30) 



The diffusion wake is a long, narrow, forward-directed stream of fluid along the negative 
axis (see flgure [T]). Its narrowness in the and directions means that, in Fourier 
space, we can take k2 and k^ large compared to ki, effectively setting k± = k. In this limit. 
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we find a simple limiting form for St- 

St-^-^(i 0)=S''. (31) 
ivki V / 

To Fourier transform S'^ to real space, we need a prescription for passing the ki = pole. 
Passing it in the upper half plane results in a stream of energy along the negative axis, 
while passing it in the lower half plane results in a stream of energy along the positive x^ 
axis — ahead of the quark. This latter possibility is unphysical in a steady-state solution. We 
find 



5f(0,f) = -^ / — .e 



iv J (27r)^ ki + ie ^32) 
'^e{-x^)6{x^)6{x^), 



where the integration is along real values, e > 0, and 

for x > 




(33) 

for X < 0. 

In [21 E], it was observed that for the case of A/" = 4 super- Yang-Mills theory, effects of order 



0{1) lead to extra terms which correspond to replacing e in (32) by ^/c^. This replacement 



broadens the wake from a singular structure along the negative xi axis to a stream inside a 
parabolic surface of revolution, as shown in figure [Tj 



The energy carried by the narrow stream described in (32) per unit time is given by 



P'^ = js^.da = j S-da = ~'^. (34) 



The first surface integral in ( 34 ) is over the boundary of a large volume enclosing the moving 



quark, and da is an outward pointing area element. This integral can be evaluated im- 



mediately, using (32), to give the result —ir^/v. In writing the second integral, we aim to 



reproduce the same result by integrating the full Poynting vector S over a carefully chosen 



surface A. From (32) we might expect that A can be any sufficiently small surface inter- 
secting the negative real axis. But in the full Poynting vector, there is viscosity broadening, 
and there are also non-hydrodynamical effects close to the quark [6l [TJ [33j. A good choice of 
A is shown in figure [2| where one takes ^1 ^ 00 with ii ^ ip ^ a/^iTJ. The large £1 limit 
insures that linearized hydrodynamics is valid, and the specified range of ip insures that the 
integration over A picks up the entire contribution from the diffusion wake and nothing else. 
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Figure 2: The diffusion wake and the integration surface A used in equation (34). The 



inviscid hmit of the diffusion wake is an infinitesmially narrow, forward directed stream of 
energy, as shown in hghter blue. After viscosity broadening, the diffusion wake thickens to 
a parabolic shape. 



The area element da points toward negative xi to agree with directionality of da in the first 



integral in (34). Combining (13), (30), and (34), we have 



(35) 



This is the result we advertised in (|2|. In section 3.3 we will describe another way of 
evaluating P*^ starting from the divergence of the stress tensor. 



3 Sound modes 

To describe the sound modes, we need to know the first column of the stress energy tensor 
in the small momentum limitj^ For purposes of simplicity, we restrict attention to a single 
scalar, so that 

i^<^ = -^G^"5^09,0-\/(0). (36) 



Referring to (20), we see that to determine the Poynting vector, it is enough to understand 
the asymptotics of E2. This field couples non-trivially both to the scalar and to Ei, E^, 
and i?4. However, in the small momentum approximation a simplification occurs: two of the 

''Near the quark, the flow is not hydrodynamical. A near-field account can probably be given, along the 
lines of 01711^. 
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equations of motion decouple, giving 



dr 



{E'^a') = 



dr 



{E'Jia^) = . 



The other equations of motion, also at order 0{l/k), are 



3a' h' 
a h 



E'i 



3a' 3h'\ ^, 



h 



'3a' ^h'\ 1 a^V'i^) 



The first order constraints couple together all the i^j's and the scalar: 



h' 



a 

Ah-E'. 
a 



2h' + 8h 



a 



a 



^^3 + -<P'l 

2h ^ 2^ 
'6'h + 2a^V'l 



^, , 3k^y\ ^, ^, i/i' k^yh' ^ 



ih' 

2~h 



kiha^ 




ikiTT^ 
k'^a^h 




0. 



V 

J 



(37) 
(38) 



(39a) 
(39b) 
(39c) 



(40a) 
(40b) 



(40c) 



where we have again written only the leading terms in a 1/k expansion. Solving (37) and 



(38) together with the requirement that the metric fluctuations vanish at the asymptotic 



AdS boundary, we get 



Eo 



E, 



62 



64 



62^— IT [h - 1) 



a(r)3 '^^hial 



dr 



Q /i(r)a(r)3 ^hia^ 



Inh . 



(41a) 
(41b) 



where we have used the equation of motion for h (7a) in the equalities on the right hand 



side, and defined hi = h'{rii). To avoid a logarithmic singularity at the event horizon, one 
must set 64 = 0. 



The alert reader will have noted that (41b) is in precise analogy to (29). The horizon 



boundary condition 64 = is likewise in analogy to ^2 = 0. The justification of these 
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boundary conditions is actually somewhat subtle. At finite k, the possible leading near- 
horizon behaviors of E4 and D2 are 

^±ivki/hr^ The minus sign corresponds to an outgoing 
wave and the plus sign corresponds to an infalling wave. Usually, the appropriate boundary 
condition is to suppress the outgoing wave. But in a small k expansion, h™'^^^^'^ = 1 + 

log h. The relative factor of ki means that the logarithmic term cannot be visible in 
a leading small k treatment. Another way to understand why one must forbid both the 
infalling and outgoing waves is that the leading order small k analysis does not know about 
viscosity, so it also doesn't know about the arrow of time; thus the boundary conditions, at 
this order, should be symmetric between infalling and outgoing waves. 

Although it is unimportant for our subsequent analysis, we note for completeness that 
the A equation of motion completely decouples from the other modes and the scalar fields 
for all values of k. At order it takes the simple form 

^(AW)=0. (42) 
dr ^ 

As with D2 and £^4, the horizon boundary condition forces A = 0{1). 



3.1 Static peturbations of the horizon 



The equations (37)-(40) are too complicated to solve in general as we did for the D set. 



Furthermore, equations (39) describe not only the response of the metric to the string, but 
also static deformations of the black hole background. Although such deformations are not 
directly related to the response to the moving quark, it is helpful to describe them explicitly 
on the dual gravity side as a warm-up to the treatment of sound modes. More precisely, we 
want to consider perturbations that are invariant under translations in the R^'^ directions 
and rotations in the dimensions. This allows us to set Hu = H22 = -f^ss and to set all 
off-diagonal components of the metric to zero. So E2 = E^^ = 0, while Ei and E^ may be 



defined as in (19). It is straightforward to show that the equations of motion take the form 



(39) and that the zero-energy constraint (coming from the rr Einstein equation) takes the 



form (40b). The expectation is that there is only one allowed deformation, and it corresponds 
to uniformly changing the temperature in the gauge theory. 

These equations are still impossible to integrate explicitly for general V{(f)). However, 
one may perform a near horizon analysis by switching the independent variable from r to h 
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and expanding 



a = ao + aih + a2h'^ + ■ ■ ■ 
r = th + Tih + r2h? + . . . 

= 00 + 4>lh + 02^^ + . . . . 



We find six linearly independent solutions: 



(43) 



A: 
B: 

C: 
D: 

E: 
F: 



El = 
Ei = - 



riaiV\<Po)h E, = -rfa^V'(0o)/i 



riaiV\<P,)h Es 
X log/i 



X log/i 



El 
El = 



1 

7^ 



E, = 



E, 



50 = 
50 = 

(50= l + r2«2^"(0o)/i 

(50= [l + r2«^^"(0o)/i] 
X log/i 

50= ^r2a2^'(0o)v^ 



-^r?«^V^'(0o)/.^/^ 



(44) 



Solutions A and B are exact, while the others, as expressed in (44), are accurate only to 



the first two non-trivial orders in small h. Applying the constraint equation (40b) to the 



solutions (44), one finds that the D solution is disallowed, and a relation is enforced between 



the E and F solutions. The upshot is that for small h, 

El = + ei - '^rlalV'{<j)o)6<j)Hh + ... 
2 

E3 = 63- -eirlalV{(f)o)y/h - 6(f)HrlalV'{(f)o)h + ... 
50 = S<Ph + '^e,rlalV'{<Po)Vh + 6<PHrlalV"{<Po)h + . . 



(45) 



where ei, ei, 63, and 50// are integration constants. In principle, three of these constants 
can be fixed in terms of the fourth by integrating Ei, E3, and 50 out to the conformal 
boundary of AdS^, and imposing boundary conditions there corresponding to the absence 
of any deformations of the lagrangian. This means that the boundary conditions on the 
asymptotic AdS boundary are such that Ei and E3 are (9(r^) near the conformal boundary, 
while 50 is 0{r^), where A is the dimension of the operator dual to 0n Once these conditions 



^There is a subtlety when A is between 2 and 3: different QFT's exist, which are conformal at least in a 
large N limit, which allow A to be exchanged with 4 — A [SJ. The simplest assumption is A > 2, but the 
precise form of the boundary condition on dcj) does not matter in our analysis. 
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from the conformal boundary are imposed, what remains is just one static deformation of 
the black hole, as expected on general groundsj^ 



For the AdS's-Schwarzschild background, the differential equations (39) can be integrated 



out to the conformal boundary of anti-de Sitter space, and the boundary conditions there 
can be imposed explicitly. They are ei = — ei = — 2e3, with no constraint on 5(f)H because 
the scalar and graviton perturbations decouple]^ To show all this, simply note that the 



terms shown explicitly in (45) provide an exact solution for AdSs-Schwarzschild. For other 



backgrounds, one usually cannot explicitly perform the integration from the horizon to the 
conformal boundary. Generically, the scalar and graviton mix, and the three constraints at 
the conformal boundary involve all of ei, ei, 63, and 5(f)H- 

Since we can not integrate the equations, we proceed by expressing the various integration 



constants in terms of the entropy and temperature. The solutions (45) imply that the metric 
fluctuations are given by 



00 



2/72, 
--eiV/i + -(es 



ei 



Hn = -63 



0{h 



1/2N 



(46) 



(recall that Hu = H22 = -f^sa-) Note that despite the l/\/h singular terms in (45), the 



fluctuations (46) are finite at the horizon. In (46) we have suppressed all terms beyond the 



ones that affect the variation in the entropy and temperature. The entropy is straightforward 
to compute: the perturbed metric to order n\ takes the form 



a 



■{h - KlHoo)dt'^ + (1 + KlHii)dx'^ + 



dr"^ 
~h 



(47) 



So the entropy per coordinate volume of x is 



1 + 4Hnf' 



(48) 



evaluated at the horizon. Thus 



5s 



7;Hu = i^^e^ 



(49) 



Determining the variation of the temperature is a bit more subtle because Gu and Grr depend 



^Intriguingly, we do not need to impose any boundary conditions at the horizon as we have in the 



discussion following (41 1. It seems that at order 1/fc the constraint equation (40b I is enough to ensure that 



Sip and all the E'i's are finite at the horizon. 

^The AdS's-Schwarzschild case is special in that a constant shift of the dilaton is an exact deformation. 
So only two additional constraints arise from the conformal boundary. Usually there are three. 
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non-analytically on h after the perturbation. One way to proceed is to change variables from 
r to f in such a way that the perturbed metric becomes 



0? 



7.2 



-hdf + dx'^ + ^ 
h 



(50) 



One easily finds 



a = a{l + -^Hu 



2 



h = h - kIHqo - nlhHii 



(51) 



dr 



4Hi 



1 ) 



dr 2 h 

all up to 0(k|) corrections which we ignore. The Hawking temperature of the unperturbed 
background is easily seen to be l/Airri. Suppose we expand f in a power series in h: 



(52) 



f = fn + fih + r2h^ + . . . . 
Then the Hawking temperature of the perturbed background is 

T- ' 



Therefore 



6T 



1 



ri 
ri 



lim 

^ h-,0 



h 



dHu 



261-63 



dh 2h dh 

According to simple thermodynamic arguments summarized in 
chemical potentials present, the speed of sound should be 



^2 ^ d\ogT ^ 1 / ei _ _^ 
dlogs 3 \e3 



(53) 



(54) 



as long as there are no 



(55) 



As a check, one may subsitute Ci = 2e3 for AdS-Schwarzschild and recover = 1/3. It will 
also prove interesting to note that the first law of thermodynamics predicts 



6S = T6s 



—-nil 



an 



-63 



4ri 2ri 

where £ is the energy density and we again assume that no chemical potential is present 



(56) 
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3.2 Small momentum asymptotics 



A striking feature of the full 0{l/k) equations (37)- (40) is that the equations of motion for 



El, E^, and 6(f), together with the zero energy constraint (40b), are the same as for static 



deformations. The two extra non homogeneous constraints (40a) and (40c) will determine 



62 and one other integration constant. Heuristically, this means that we can describe sound 
waves as a modulation of the static deformation, plus some profile for E2 that incorporates 
non-zero Poynting vector. Technically, it means we can set 



ei 



;i + 3ae3. 



(57) 



This is not really a constraint on ci , instead we are simply trading ei / 63 for = d log T/d log s 
as a parameter for describing temperature-deformations. Only by fully integrating the dif- 



ferential equations (39) and (45) and imposing appropriate constraints at the conformal 
boundary could we determine c^. Plugging (45) and (41) into (40), and imposing 64 = 



as well as (57), one obtains a pair of linear equations in 62 and 63. Solving these equations 
leads to 



62 



iki v'^kf — c^/c^ 



63 



2riki'K^ 1 



(58) 



Using (30) and (58) in (22) we obtain the Poynting vector to leading order in small k: 



Si 
Si 



icIkiTT^ 1 + v"^ 



ivki 
icj.k±n^ 



V v'^kf 
1 + 



j'^kf — c^k"^ 



c^k"^ 
+ 0{l). 



+ 0{l) 



(59) 



When Cs = l/V^, (59) reduces to the leading order results found in [2]. It is natural to 
identify the contribution of sound modes as 



S' = (^Sl SI 



where 



SI 



h 
ki 



S' 



Defining S as in (31), we see that S 



V 

S' + S'^. 



j'^kf — c^k"^ 



(60) 



(61) 



In order to completely characterize the energy dissipation, we need also to compute 
the fluctuations in the energy density. This is difficult to work out in a general setting 
because the Brown- York tensor probably does not correctly capture (Too) owing to effects 
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from holographic renormahzation. Also one apparently has to solve for Ei and out to 
the boundary, which seems impossible. But there is a simple shortcut — though admittedly 
heuristic: it is simply to use the thermodynamic relation 6S = T6s in the k-th Fourier mode. 



This should work at leading order in small k. The result, carried over directly from (56), is 



6S 



'2ri 



eg = ikiTf^ 



v'^kf — c^k"^ 



(62) 



When Cg = l/V^, (62) reduces to the leading order results found in [2]. 



3.3 Non-conservation of the stress tensor 



The results (59) and (62) are enough information to determine the first row of the stress 



tensor completely at small k. More precisely, consider the subtracted stress energy tensor 
expectation value at time t = (still in the rest frame of the plasma): 



(T""(0,f)) - (T'"") 



bath ; 



(63) 



where (T'""')bath is from the infinite, static, thermal medium. What we have computed, in 
Fourier space, with k^ = and k2 = kj_ > 0, is 



gr^Om ^ 1^^^ 



(64) 



(Recall that we use S for the response of the Poynting vector to the string instead of SS.) 
Rotational invariance around the 1 axis determines ^T^"* for other values of k. As discussed 
in [2], |32], the failure of the conservation law dmST"^"' = is a good measure of energy 
dissipation from the quark. The non-conservation law can be expressed in Fourier space as 



ikm.ST 



Om 



pt 



(65) 



where 



(^—vki ki 



(66) 



(setting ko = —vki comes from the co-moving ansatz.) The right hand side of (65) was 



computed in (13). It is the total rate at which the quark deposits energy into the plasma: 



a positive quantity. A check on our computations so far is to plug (59), (62), and (64) into 



(65) and check that the right hand side agrees with (13). It does. 
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As in P], we associate 5£ entirely with dissipation into sound waves. This makes sense 
because, in a hnearized hydrodynamical analysis, the diffusion wake shows up only in the 
Poynting vector. Thus we have a splitting 



where 

5t!Z.a=U£ si si o) 

(68) 

^^dZse = (o sf o) . 

There is an associated splitting of the non-conservation effect: 

•T rpOm _ ps 
'"-m-t sound ~ ^ 



■7 rpOm _ pd 
t^m-i diffuse ~ 



(69) 



We regard (69) as a definition of and P'^. It is immediate from the preceding discussion 
that 

P"" : P'^ : P' = 1 + : -1 : v\ (70) 
which is the result ([s]) that we advertised from the start. 



It is clear that in splitting up ST as in (68) we are making choices rather than doing 



something inevitable. To justify this splitting we note that in linearized hydrodynamics, 
neglecting viscosity, the Fourier space form describing sound waves involves u"^ — c^k^, while 
the analogous form for the diffusion wake is just u (see for example [l]). This precisely 



coincides with our choice of sound and diffusion modes in equations (31 ), (61 ) and (62). The 



diffusion part of the stress tensor is exactly the long, narrow stream of plasma behind the 



quark, which we justified in a different way in the discussion leading up to (31). 

Furthermore, when v > Cg, the split we use leads to exponential suppression of all com- 
ponents of ST^^^^ and ST^^^^^ for > 0, which is sensible because this region is causally 
inaccessible to sound waves emanating from the quark. Another natural split is to associate 
contributions to the stress tensor from the D set purely with the diffusion wake, and contribu- 
tions from the E set purely with sound modes. This is the same transverse/longitudinal split 
discussed in [1| and it leads to power law tails of components of the stress tensor cancelhng 
delicately between diffusion and sound modes. 

It is interesting to note that when 0, the total perturbation in the stress tensor 
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becomes 



vki 



1 



with our usual convention that = and k2 = k± > 0. In position space, 



1 



(71) 



(72) 



which describes energy deposition and flow only in the path of the quark. There is of 
course some viscosity broadening, presumably visible at higher orders in k, but no pressure 
broadening — because there is no pressure. Sound waves don't really exist in this limit. But 



the diffusion wake, as defined in (31), is completely insensitive to the speed of sound. Its 



form and its strength appear to be a universal feature of string theory constructions starting 
from Einstein gravity coupled to strings and scalars. 



4 Relation to linearized hydrodynamics 

The stress energy tensor of the thermal plasma is conserved except at the location of the 
quark: 

d^ST"''' = r5{x^ - vt)5{x^)5{x^) . (73) 



In Fourier space. 



where 



ik^ST'^'' = J" 



km = (^—vki ki k± 



(74) 



(75) 



as in (66). The n = component of (74) is precisely (65). From (12) and (13) we read off 



J" = [vn^ TT^ o) . 



(76) 



We emphasize that our derivation of (74), for n = 0, did not directly use a hydrodynamic 



approximation: it came from the small k limit of an AdS / CFT computation. The purpose of 
this section is to understand how the AdS / CFT computation matches with linearized hydro, 
including (for A/" = 4) a sub-leading effect in small k. 

The constitutive relations for linearized hydro amount to 



5T, 



hydro,ij 



^s^^ij ^kiSj -\- ikjSi ^^ijikiS^ , 



(77) 
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where 



"^-^ hydro 



>Ji — hydro ■ 



(7J 



(Recall that we use mostly plus signature.) The equations of motion for linearized hydrody- 
namics are 



i Aim, hydro 



jn 

hydro ' 



(79) 



But now there is no reason why the source term J^ydro should be pointlike. Instead, in 
position space, J^ydro probably has width comparable to Tg (the attenuation length) and 
must be parametrized in some fashion. The authors of [H ES] used an equivalent of 



>^hydro = (eo 9o + ki9i k±gi . 



(80) 



In [H EB] two scenarios where considered. In scenario 1, Cq = vg^ is a Gaussian in real space 
with width F^, while gi = 0. Scenario 2 corresponds to the opposite extreme, Cq = go = 
and gi a Gaussian with width Tg. 

A natural though non-unique prescription in matching our results to linearized hydro is 
to set 

ST^l, = 5T^^ (81) 
for m = 0, 1, 2, 3 and then to use (77) and (78) to determine the other components of ^T^^roj^ 



If we follow this procedure to order 0{k ^), see (59) and (62), we find 



irrpmn 
f^-' hydro 



1 _ ic^ki(l+v^) 
ikiv D(cf fc2— fcjti^) 

_ ic^kjjl+ip)_ 
v{c'^k'^ — k^v^) 





1 _ ic^Mi'^+f) 
ikiv v(c^k'^ — k'fv'^) 
ik-i ci 



-c^k'^+k'fv'^ 






tc^fcx(l+t;^) 
v{c'^k^ — k'fv^) 



ifeic2(l+i;2) 











— c^k^+kjv^ / 



(82) 



If we plug (|82|) into (|79|), we obtain an expression for J^y^ro coinciding with (|80|) with 



eo = VTT^ + 0{k) (70 = TT^ + 0{k) gi = 0(1) . 



Fourier transforming the leading terms shown in (83) leads to position results which are 



delta-functions, as in (73). Because of the explicit factors of in (80), non- vanishing eo and 

^One possible alternative prescription for determining ST^'^^^ is to require that at order 0{k), J^ydro 
only longitudinal components. We've checked that the results we describe below remain unchanged under 
such a redefinition. 
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go tend to dominate over non-vanishing gi at sufficiently long distances. More specifically, 
given Co, go and gi in Fourier space, one may define a dimensionless figure of merit. 



7i 



1 
R 



9i{k = l/R) 



eo{k = l/R) 



(84) 



such that 7i <^ 1 favors scenario 1 and 71 ^ 1 favors scenario 2. i? is a typical distance 
at which one is interested in the strength of the fiow induced by Cq and go as compared to 
the fiow induced by gi^ (We assume that cq and go have similar magnitudes.) If Jhydro 
localized at a scale Tg, then for a source such as the one we have described, where cq and go 
are non-vanishing at leading order, one expects 



71 



R 



(85) 



Because one only trusts linearized hydro for R ^ Tg, our results correspond to scenario 1. 
To get scenario 2, one needs 71 ^ 1. 



The conclusion (85) was based on dimensional reasoning and so should be treated with 



some caution. We can calculate a precise value for 71, as well as more information about the 
structure of the source terms Cq and go, by calculating the stress tensor at the first subleading 
order in k. That is, we expand 



5T' 



jrrjimn 
"-'hydro 



+ 5r(i)™" + C(A;) 



ST 



(0) mn 



(86) 



+ 5Til\Z + Oik) 



and impose (77) and (81) at order 0{1) to get ST^^^^"". These subleading contributions 
to the energy momentum tensor are currently available only for the J\f = 4 theory [32j. 
Adapting the notation in [32] to the current conventions, we find that 



ST^') = STi^lZ + d^^Sona^{o -2 1 1 



3^' 



^7) 



This corresponds to 



J"" = (v 1 - 2vkiTs vk±Ts o) vr^ + 0{k^) , 



^J. Casalderrey-Solana has suggested to us that 71 could be changed by sending R VT^R in (84 1. This 



results in a rough expectation 71 ~ y^Vj/R in place of ( 85 ) 
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where we have used = l/SvrT. From (88) we read off 



eo{k) = TT^v + 0{e) go{k) =n^{l-3vhTs) + 0{k'') gi{k) = n/rvTs + 0{k) . (89) 

These results agree with the expectation that Jhydro structure on the scale F^, and it 
implies that 

71 = ^ + 0{R-') . (90) 

Taking F^ = 0.07 fm and i? = 7fm to make a rough comparison with heavy ion collisions, 
one finds 71 = 0.01. 

It seems likely that 71 will be be small for phenomenologically interesting values of R even 
for AdS/CFT constructions including scalars, unless perhaps the scalar lagrangian has some 
very special or extreme form. Nevertheless, it would be interesting to see what form these 
subleading expressions take in, for example, the Af = 2* theory [23 ESI EI], or the casacading 
gauge theory [3E1 EE], EH EH] . It seems to us more likely that 71 could be altered significantly 
by considering a different source from the heavy quark we have employed throughout. 



5 Jet-splitting and the challenge of diffusion wakes 

In this section we attempt to bridge between our results from AdS / CFT and the phenomeno- 
logical literature on heavy ion collisions. In doing so we should bear in mind all the usual 
caveats to such comparisons. In particular: the AdS/CFT results apply in the limit of large 

and large 't Hooft coupling gy^N; the quark in our treatment is infinitely heavy and has 
constant velocity; and the holographic dual to QCD is not known. 

In heavy ion collisions, a well studied experimental probe of the interaction of a hard 
parton with the medium is the relative azimuthal angle between two energetic hadrons 
emitted from the interaction region. Histograms of hadron pairs always show a peak at 
A(f) = 0, most likely meaning that the two hadrons were part of the same jet. In proton- 
proton collisions and sufficiently peripheral heavy ion collisions where the extent of the 
medium is small, another peak at A0 = vr arises, owing presumably to events where two hard 
partons are created with back-to-back momenta — at least, back-to-back in the azimuthal 
direction. See figure [3j 

Jet-broadening or jet-splitting refers to the change in shape of the peak at A0 = vr. Di- 
hadron histograms published by the PHENIX collaboration in [39] show a definite double- 
peaked structure: a minimum at A0 = tt and peaks at A0 ^ vr ± 1.2. On the other hand. 
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Figure 3: A cartoon of a hard process in a heavy ion colhsion. The blue and green crescents 
are the parts of the nuclei that didn't interact. The gold region is the thermalizing medium. 
The red dot is the vertex of a hard process that occurs early in the collision, producing 
two partons with back-to-back momenta in the azimuthal direction. For simplicity we show 
all particles in the x-y plane meaning at mid-rapidity. If particles 1 and 2 escape to form 
hadrons, it would correspond to A0 n. If instead particle 2 stops in the medium but emits 
a high-angle secondary which is then observed as 3, then is significantly different from 

TT. 
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the di-hadron histograms pubhshed by the STAR collaboration in [10] do not exhibit a 
pronounced double-peaked structure and are discussed in that paper in terms of softening in 
PTi broadening in A0 and rapidity, and thermalization. Momentum cuts were significantly 
different in these two studies: 

STAR : 0.15 GeV < < 4 GeV < p*"^ < 6 GeV 

(91) 

PHENIX : 1 GeV < p^"*^"' < 2.5 GeV < p^"^ < 4 GeV . 

star's pseudo-rapidity acceptance, \rj\ < 1, is also significantly greater than PHENIX's, 
|?7| < 0.35. Subsequent experimental analyses, including [HI HJl |l3l EH US], support the 
conclusion that at least for intermediate momentum hadrons, similar to the PHENIX cuts 



shown in (91), splitting of the away side peak does occur. 

A candidate explanation for jet-splitting is that a hard parton traveling through a quark- 
gluon plasma (QGP) loses energy in large part by producing a sonic boom [T],|l6]. The sonic 
boom would then propagate through the plasma until freeze-out converts it into an excess 
of hadrons emitted in roughly the direction of the Mach angle, relative to the direction of 
the original parton. There are theoretical alternatives (see for example jlZl HHl HHl [50]), but 
our string theory results bear less directly upon them. In [Ij, two scenarios were presented 
for the distribution of energy into sound modes and diffusion modes. In section |4] we have 
shown that in scenario 1 the relative strength of the sonic boom and the diffusion wake is 
the same as in our string theory analysis. After Cooper-Frye hadronization on a surface of 
fixed time, it shows no jet-splitting due to the relatively strong diffusion wake. In scenario 2, 
only longitudinal modes are sourced by the energy momentum tensor so the diffusion wake 
is entirely suppressed and a jet-splitting effect is recovered after Cooper-Frye hadronization. 
Further work [36] showed that for scenario 2 to match to data, the energy loss per distance 
needs to be chosen quite large — perhaps unrealistically so. Another study [51] included 
expansion of the plasma and used source terms similar to scenario 1. It also concluded 
that a minimum at 6(f) = it is hard to achieve. In [521 ES], a model including Mach cone 
propagation through the QGP has been argued to fit experimental data provided the fraction 
of energy going into the sonic boom is 0.75 or more. But it is not yet clear to us how this 
fraction is related to the relative strength of the sonic boom and the diffusion wake as we 
compute it. 

The robustness of the diffusion wake in the string theory construction we have explored 
favors scenario 1. It is therefore a concern that this scenario seems less likely to agree with 
data. It is natural for us to hope that a diffusion wake really is generated when a hard parton 
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traverses a quark-gluon plasma, but that it is erased or disguised in the data. More reahstic 
theoretical treatments should include expansion of the plasma, as in [5ll [521 [53] and perhaps 
also a more sophisticated treatment of hadronization; as compared to the sonic boom, the 
diffusion wake is deeper inside the expanding medium, so it has more chance to broaden, 
soften, and thermalize before freeze-out. Some hints of this were observed in pT]. There 
also issues of detector acceptance which may be relevant, see for example the discussion in 
[53] . Finally, we should keep firmly in mind our starting assumption that the quark is very 
massive. If the trigger particle could be tagged as originating from a heavy quark, then the 
distribution of associated hadrons would be more directly related to our analysis. Studies 
cited thus far do not include flavor tagging; see however [5^ . 

In summary: We have argued that the strength of the diffusion wake relative to the 
drag force is universal in AdS/CFT constructions based on a trailing string in a background 
supported by arbitrary scalar matter. This universality is in some ways comparable to the 
universality of the rj/s ratio |55l [561 EI] (see also [58]): it arises because of a decoupling 
of certain metric perturbations from the scalar dynamics, and it may suffer corrections in 
inverse powers of N and Qym^ ■ Also, it is worth remembering that our slightly weaker 
result on the universality of dissipation through sound modes does not specify how well 
this dissipation is focused at the Mach angle. Phenomenologically, it may be challenging to 
accommodate as strong a diffusion wake as we predict. 
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A Holographic renormalization 

The prescription for calculating the expectation value of a gauge theory operator dual 
to a field $ in the supergravity limit of string theory is to identify the generating function 
of with the supergravity partition function [321 [HO]- Sources for the operator 
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in the gauge theory are identified with boundary values of $ (modulo some issues with low 
dimension operators [34].) In the saddle point approximation, calculating one point functions 
reduces to varying the supergravity action with respect to the boundary values of the fields 
$ near the AdS boundary. 

<o*>=lL"J^a|5ySscc«.. (92) 

Here g is defined through ds'^ = ^ (dr^ + gmndx"^dx"'), and g^^^ is the leading 0{r'^) term in 
a series expansion of g. Similarly, ^'^^^ is the leading 0{r^~^) term (with A the dimension 
of 0$) in a series expansion of $. 



At this point the prescription (92) is not well defined since generically (92) will give di- 
vergent results. Simply ignoring these divergences may lead to wrong Ward identities. One 
prescription for systematically dealing with these singularities is holographic renormalization 
[6T| [62] (see for a review). This prescription seems to work well even for non asymp- 
totically AdS backgrounds [64]. It comprises of adding a boundary action to Ssugra- 
Set may be composed only of the boundary values of the field $, and is constructed so as to 



precisely cancel the r ^ divergences which appear in (92). 

In the current setup, we are interested in the expectation value of the energy momentum 
tensor 

2 d 

(Trnn) = hm — — {SsUGRA + S^t) (93) 

where is the Brown York stress tensor, which is what we would obtain by varying the 
on-shell Gibbons Hawking term in the action and simply discarding the singularities. More 
precisely, we are interested in the contribution of the trailing string, or moving quark, to the 
energy momentum tensor to linear order in k^, i.e. 

((5T„„) = 5T^J + 5X„„. (95) 

The purpose of this section is to show that SX^n is diagonal and to evaluate the off diagonal 
parts of (STmn) in terms of the near boundary value of the metric fiuctuations. 
Expanding the various scalar fields in a power series in r, we find 

0^(r) = r^-^^(0^(o) + rV'^'^ + •.•)+ r^^(0'(°) + . . .) (96) 
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where A/ is the dimension of the operator dual to 0^. If Aj is an integer then additional log- 
arithmic terms will appear in the series expansion for (p^ . Fields dual to irrelevant operators 
have A/ > 4. In that case the power series expansion in r is divergent at the boundary and 
the sources 0^ should be taken to be infinitesimal, so that they vanish on-shell |60l [61] . 
For clarity, we shall omit the index J, and reinsert it only when required. The coefficients 
0^"-' in (96) can be related to (p^^^ by solving the equation of motion for order by order in 



r. The coefficients 0'-"'' are related to 0*^°^ in a similar manner. 0'^°-' and 0*-'^'' are determined 
through the boundary conditions. For the case at hand the response of the field to the string 
is such that the leading near boundary contributions vanish. Using the same notation as in 



(14) 



^ + 50 (97) 
we have (50(°) = 0. Also, in our construction the background value of the fields depend 



only on the radial coordinate, meaning that in the notation of (97) 0'-^ is a constant. A 



similar analysis follows for a near boundary expansion of the metric 

^ ^(0) , (2) 2 , -(0) 4 , /QgN 

Consider now the boundary counterterm action Set- As discussed earlier it must be com- 
posed of the near boundary metric 7mn(e) = '^fi'mn(e) and fields 0(e), where e is the radial 
coordinate which will be taken to zero at the end of the calculation. Writing the metric 
dependence explicitly, we find 

Sct = ^ ! VlLctd'x = -i^ / V7 (^A + RLb + Y'L.j + 0{e')) d'x, (99) 

where the integrals are over a surface which is a distance r = e from the asymptotically AdS 



boundary. In (99) we have used 



7mn = e-' 77-6-^ ^™"^e2 i?™„fc;~e° (100) 

i?mn~e° ^2~e' D-e' (101) 

and we've also rewritten any n(f)^Lj expressions as 'y"^"'dn(p^dmLj. 
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We wish to evaluate the finite terms in 



2d L"^ 1 d 

- {Set + Sb) = lim -2-^ — {Set + Sb) 



(0) c}^{0)m" 



'9 



e^o nle^ ^7 d-f 
L2 



hm 



2^2 



'~1mn 

{L,t + Lb) + 2 



dL 



ct 



(102) 
(103) 



where Sb = J ^J^Lbd^x is the on-shell boundary action which gets contributions from the 
various scalar fields. 



We start from the last expression in the parenthesis in (103). Using (99), the leading 
contributions to which come from the trailing string will look like 



e ^ {5{RmnLR) + 5Lmn + •••)• 



(104) 



The second term in (104) vanishes since it must involve two derivatives in a direction tangent 
to the radial AdS direction. One of these will always act on a background field (p^^^ which 
depends only on the radial direction, giving dicj)^^'' = 0. 



Consider the first term in (104). and 6Lji are polynomials of the scalar fields with 
LptjSLji ~ e" with n > 0. (The inequality for n follows by considering (96). If n < then 
Lri^Lr ~ e" would imply a non vanishing source term for a non-normalizable operator.) 



So we can expand Lr = ^s=o-^R^^'^ ^^'^ ~ ^s=o^-^r'^^- Similarly Rij = for i 7^ j 



because the background metric g is independent of the Xj directions and 6Rmn = Ss=4 
since the boundary conditions on the response of the metric to the string are that the metric 
vanishes close to the boundary]^ This brings us to the conclusion that 6{RijLji) ~ e'^ so the 



first term in (104) will not contribute to the stress energy tensor through (99) and (103) 



Finally we consider contributions coming from j^nn [Lb + Ld)- Expanding 



In 



(g(0) +0(2)^2 -(0) 4^ N 
\^mn ' Hmn^ ' Hmn^ ' / 



(2) ,2 , -(0) 4 



and Let + Lb = L^^^ + L^^^ + O(e^), we see that any finite non diagonal contributions of 
^7jj {Let + Lb) to {Tmn) can only come from (jmn (^L^a + -^1°''^ But such expressions will 
also contribute to a divergent C(e~^) term in ^7 (La + Lb) which will not be compensated by 



^In |32j it was shown that even when the boundary conditions imply that the metric fluctuations vanish 
at the AdS boundary, the stringy source induces an 0{r^) k independent term in the series expansion for 
the metric fluctuations. This 0{r^) contribution is apparently related to the infinite mass of the quark and 
corresponds to a delta function at the location of the quark when going to real space. Regardless, it does 
not appear in the 1/fc expansion that we are considering here. 
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the Gibbons Hawking term since it involves the scalar fields. Therefore such a contribution 
will not exist. 

We are left with evaluating the Brown York stress tensor. In a general setting, this tensor 
is defined on a hypersurface S whose outward pointing normal is A^"^. (Outward pointing 
means pointing toward the conformal boundary.) The definition is 

T^l = Kmn - + 1) aln , (105) 

where is the induced metric on S and Kmn is its extrinsic curvature. This expression 
comes from varying the Gibbons Hawking term in terms of the boundary metric 7^™- The 
indices m and n run from to 3 (along S) while /i runs over all five dimensions of the 
asymptotically AdS spacetime, whose radius of curvature near the boundary is L. In general, 

9mn = 9mn " N^Nn Kmn = -gJ^kNn , (106) 

but in axial gauge, when S is a surface at fixed r, one has the simpler expressions 

9mn ~ 9mn K^n — / „ drQmn ■ (107) 

The metric of the boundary theory and the holographic stress tensor must be computed in 
a coordinated way, because both are affected by conformal transformations. A consistent 
prescription is 

^mn = -^gmn {Tmn) = ^ hm tt^T^J . (108) 

For the ansatz ([6|, one finds g^^^ = i] = diag{— 1, 1, 1, 1} and 

{5Tmn) = lim —{Hmn ' v'^Hi^Vmn) , (109) 
r— ►0 L 



where we have used (15). 
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